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ABSTRACT 

The  radiated  noise  of  ships  or  other  underwater  sources  are 
typically  characterized  in  terms  of  a  far-field,  plane-wave 
equivalent  source  level  ba.sed  on  measurements  assumed  to 
have  been  acquired  in  a  free-field  environment  such  as  a  deep 
water  test  range.  Measurement  of  ship  noise  in  a  harbor 
environment,  where  multiple  reflections,  high  background 
noise  and  short  propagation  paths  are  the  nomi,  violates  the 
conditions  tliat  as.sume  the  ship  is  a  radiating  simple  source. 
Careful  analysis  is  required  to  arrive  at  a  valid  estimate  of  far- 
field,  free-field  source  levels  from  such  measurements. 

This  work  presents  results  from  a  test  conducted  at  the  US 
Navy  Acoustic  Research  Detachment  in  Bayview,  Idaho  during 
summers  2010  and  2011.  A  line  of  omnidirectional 
hydrophones  was  deployed  from  a  barge  adjacent  to  a  moored 
test  vessel  to  obtain  radiated  noise  measurements  fi-om  several 
shipboard  .sources.  A  series  of  lest  signals  was  also  transmitted 
through  calibrated  acoustic  sources  to  evaluate  the  effectiveness 
of  post-processing  techniques,  as  well  as  line  array 
beamforming,  in  minimizing  reflected  path  contributions  and 
improving  signal-io-noise  ratio.  Methods  of  estimating  far- 
field,  free-field  equivalent  source  levels  based  on  these 
measurements  are  presented. 


INTRODUCTION 

Ship  radiation,  as  well  as  other  underwater  sources,  is 
typically  characterized  in  a  deep  water  environment  using  a  far- 
field,  free-field  assumption,  which  is  almost  exclusively  the 
only  method  covered  in  array  signal  processing  literature.  This 
assumption  states  that  the  ship  is  a  radiating  simple  .source  and 
the  wavefronts  received  by  the  array  can  be  considered  to  be 
planar.  Literature  typically  defines  the  far-field  .starling  at 


where  r  is  the  distance  from  tlie  center  of  the  array  to  the 
source,  L  is  the  largest  dimension  of  the  array,  and  X  is  the 
operating  wavelength. 

Measurement  in  a  harbor  environment,  however,  violates 
these  assumptions  for  frequencies  in  the  range  of  interest  and 


other  signal  processing  metliods  must  be  explored.  Therefore, 
we  consider  beamforming  to  reduce  the  multipath  contribution 
to  the  signal  as  well  as  other  background  noise.  This  also 
allows  us  to  take  into  account  the  spherical  geometry  of  actual 
wavefronts. 
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FIGURE  1.  Research  site  geometry. 


Lake  Pend  Oreille  in  Bayview,  Idaho  is  one  of  the  deepest 
lakes  in  the  country,  with  a  maximum  depth  of  approximately 
1200  feet.  It  is  also  home  to  a  US  Navy  Acoustic  Research 
Detachment  (ARD).  It  is  an  ideal  location  for  acoustic 
experimentation  due  to  small  amounts  of  boat  traffic,  making 
the  ambient  levels  quite  low.  An  array  of  14  omnidirectional 
hydrophones  was  deployed  off  a  barge  next  to  a  moored  quarter 
scale  DDGIOOO,  the  SeaJet,  in  a  shallow  (<50ft  deptli)  harbor. 
Two  sources  were  tlien  lowered  into  the  water  column  off  the 
side  of  the  barge  and  transmitted  various  test  signals.  The 
geometry  of  the  environment  is  shown  in  Figure  1.  These 
signals  and  measurements  were  used  to  obtain  calibrated  data  in 
order  to  evaluate  several  post-processing  techniques.  This 
paper  explores  various  meUiods  of  post-processing  and  their 
effectiveness  in  minimizing  reflected  path  contributions  as  well 
as  improving  signal-lo-noise  ratio. 


OMNIDIRECTIONAL  PROCESSING 

An  array  of  14  omnidirectional  hydrophones  witli  16  inch 
spacing  between  elements  was  deployed  in  both  a  harbor  test  as 
well  as  a  deep  water  test  in  tlie  lake.  Basic  signal  processing 
was  used  to  find  third  octave  band  levels.  Using  kno\\Ti  depths 
and  ranges,  expected  transmission  loss  can  be  calculated  from 
tliese  levels  using  tlic  standard  convention: 

TL  =  20  log(r) 

Signals  transmitted  in  the  deep  water  tests  hold  very  close 
to  tin’s  convention,  with  test  values  falling  within  1  dB  re  IpPa 
of  the  expected  value.  This  reveals  that  we  are  operating  in  tlie 
frcc-ficld.  with  very  little  to  no  contribution  from  reflections. 
This  can  be  seen  with  a  broadband  noise  signal  plotted  against 
the  ambient  in  Figure  2a.  Signals  transmitted  in  the  harbor 
tests,  however,  showed  transmission  loss  lower  than  tlic 
convention,  witJi  test  values  being  over  4  dB  re  IpPa  higher 
than  tlie  expected  value.  Tests  done  in  the  harbor  environment 
arc  greatly  contaminated  by  multiple  reflections  from  not  only 
the  surface  and  the  bottom,  but  also  other  underwater 
interfaces.  Tliis  can  cau.se  many  problems  in  determining  an 
equivalent  source  level.  This  is  shown  in  Figure  2b,  where  tlie 
same  source  was  radiating  the  same  power  level  over  a  slightly 
difl'erent  range.  A  very  different  third  octave  band  level,  and 
tlierefore  source  level,  can  be  observed  between  the  far-field 
measurement  and  the  measurement  in  tlie  multipath 
environment. 


Figure  2.  Broadband  noise  generated  by  a  J9  omnidirectional 
source  shown  In  one  third  octave  band  levels.  Reference 
hydrophone  one  meter  from  the  source  (green),  bottommost 


array  hydrophone  (red),  ambient  level  (blue).  The  band  of 
interest  here  is  the  4  kHz  band.  Figure  (a)  depicts 
measurement  In  a  deep-water  environment  with  the 
expected  30  dB  of  transmission  loss  while  Figure  (b)  depicts 
measurement  in  a  shallow-water  environment  with  10  dB  of 
transmission  loss  where  18  dB  Is  expected. 

LINEAR  FAR-FIELD  BEAMFORM  PROCESSING 

In  far-field  bcamforming,  we  begin  wiili  a  bcamforming 
equation  for  array  processing  given  in  [1]: 

^  A 
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From  this  base  equation,  we  make  the  a.ssumptions  lliat  all 
rt  are  approximately  parallel,  so 

r»(^N-  l)d 

and 
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for  all  /  where  r  is  the  distance  from  array  center  to  source.  We 
also  add  a  phase  delay: 

2n  2n 

(pi  =  i  -T- Ar  =  i  —d  sin  Oq  (2) 
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where  d  is  tlie  element  separation  and  Oo  is  the  incident  angle. 
With  tliese  a.ssumptions,  we  arrive  at  a  corrected  equation: 

N 

p{r,  e,  t.  0o)  =  ^  (3) 

where  Ar  =  rising. 

*niis  is  very  similar  to  equation  (7.8.2)  in  [1]  except  with 
an  added  phase  delay.  Strict  time-domain  bcamforming  is  not 
an  option  with  this  set  of  data  because  our  resolution  is  not  fine 
enough  to  establish  a  proper  angle  of  incidence.  So  instead  w  e 
turn  to  frequency  domain  bcamforming.  A  Fourier  Transform 
on  tin’s  Eqn.  (3)  gives  us  the  equation: 

This  is  only  multiplied  by  tlie  unsteered  Fast 

Fourier  Transform.  If  we  look  at  equations  given  by  [2],  we 
can  write  this  as  a  vector  equation: 

Yif)  =  d^'(ej)Wx(J) 


(5) 


Here,  0  is  the  incident  angle,  W^is  a  weighting  matrix: 


/)  is  a  steering  vector  which  depends  on  incident  angle, 
frequency,  element  spacing,  and  speed  of  sound: 


and  .v(/)  is  the  data  vector: 


We  can  use  this  knowledge  to  apply  a  far-field  beamforming 
algorithm  to  our  data,  lliis  allows  us  to  create  a  bearing-time 
record  of  a  specific  run  over  a  certain  period  of  time. 


Nt»r*Fteld  Be«nn9-Iime  Record  of  MOo  62^ 


Figure  3.  Bearing-Time  Record  of  broadband  noise  generated 
by  a  J9  omnidirectionai  source  in  shallow  water  over  a  5 
second  window. 


Averaged  BTR  over  all  times  fer  MOo  62g  with  center  frequer^cy  3720  Hz 


Figure  4.  Average  of  Figure  (3)  over  the  5  second  time 
window. 

In  Figure  (4)  we  can  see  three  distinct  paths,  likely  the 
direct  path,  surface  reflection,  and  bottom  reflection.  However, 
we  can  also  see  that  this  is  muddied  with  other  reflections  as 
well.  We  also  now  know  the  direction  of  the  noise  source,  but 
not  the  actual  location,  specifically  the  range.  This  is  crucial  in 
determining  the  transmission  loss,  and  ultimately  the  source 
noise  level. 


NEAR-FIELD  BEAMFORM  PROCESSING 

Far-field  bcamforming,  while  quite  convenient  and 
accurate  in  far-field,  free- field  measurements,  loses  its 
precision  in  multipath,  near-field  environments  where  far-field 
assumptions  cannot  be  made.  Once  again,  starting  with 
Equation  (1),  we  can  use  some  basic  geometry  to  eliminate  the 
far-field  assumptions  and  account  for  the  curvature  of  the  wave 
front.  We  will  also  measure  from  the  geometric  center  of  the 
array,  instead  of  the  top-most  element. 


Figure  (3)  shows  us  the  strongest  angles  of  incidence  from 
the  source  (presumably  the  direct  path,  and  first  reflections  off 
the  surface  and  bottom  interface).  Averaging  these  levels  over 
the  time  window,  we  are  able  to  see  the  contributions  from  each 
path  at  each  angle. 


Figure  5.  Test  site  geometry  that  includes  near-fteid  paths. 


Figure  (5)  shows  the  geometry  of  this  near-ficld  problem  i 

and  gives  us  a  relatively  easy  method  of  determining  a  phase 
shift.  Recalling  tliat  the  phase  shift:  ^ 

— 

Tliis  lime,  wc  cannot  assume  that  all  r/  are  equal,  so  * 

where  r,  is  the  distance  from  the  source  to  each 
individual  element  and  r  is  tlie  distance  from  the  center  of  tlie  i 

array  to  tlie  source.  The  distance  /'  can  be  represented  by: 
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We  can  then  use  the  distance  from  each  individual 
element  to  tlie  center  of  the  array: 


(6) 


where  N  is  the  total  number  of  elements,  i  is  the  element 
number,  and  d  is  the  separation  distance  between  each  element, 
and  the  law  of  cosines  to  find  the  distance  from  each  element  to 
tlie  source: 


(7) 


The  angle  steering  angle  for  each  individual  element  is 
also  no  longer  constant,  and  can  be  determined  by: 


(8) 


This  can  now  be  used  in  Equation  (1)  to  find  tlie  specific 
pressure  without  making  any  far-field  assumptions.  Once  again 
wc  can  take  an  FFl'  of  this  equation  to  find: 


(9) 


Once  again,  it  appears  tlie  only  addition  to  tlie  unsteered 
FIT'  is  a  phase  delay  term,  ,  that  is  dependent  on  both 

steering  angle,  </>„  and  range,  R. 
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Figure  6.  Illustration  of  the  array  with  the  source 
normalized  to  the  origin  and  highlighted  in  green.  The  top 
hydrophone  is  highlighted  in  blue  because  it  is  shallower  than 
the  source.  This  array  shape  agrees  with  visual  observations 
of  the  array  out  of  the  water.  The  distortion  in  the  shape  of 
the  array  is  likely  caused  by  the  weight  of  the  cables  and 
hydrophones  causing  stress  on  the  array.  Units  on  each  axis 
are  meters. 

A  complication  in  this  method  lies  in  tlie  error  of  precise 
knowledge  of  the  array  element  locations.  Using  other  methods 
of  source  localization,  wc  determined  that  our  array  was  not,  in 
fact,  perfectly  linear.  Instead,  it  more  resembled  tlie  shape 
shown  in  Figure  (6).  This  shows  a  max  difference  of  nearly 
one  meter  between  the  bottom  element  and  one  of  the  topmost 
elements.  With  wavelengths  in  our  range  of  interest,  this 
location  difference  can  change  the  initial  phase  of  the  signal 
enough  that  it  voids  the  phase  shift  that  we  are  trying  to  apply. 

Tlicse  have  not  yet  been  implemented,  however,  much 
research  suggests  that  near- field  beam  forming  will  reduce  side- 
lobe  contribution  even  more  and  allow  us  to  determine  the 
source  noise  level  with  the  most  accuracy. 

CONCLUSIONS 

Underwater  noise  source  radiation  is  difficult  to 
characterize  in  a  shallow  water  environment  where  multiple 
paths  contribute  largely  to  the  overall  noise  level.  We  have 
presented  various  degrees  of  measurement  and  processing  skills 
that  ultimately  may  be  able  to  accurately  determine  tlie  noise 
source  level  of  an  object  in  tliis  environment  when  notliing  is 
known  about  tlie  object.  Further  investigation  will  determine 
just  how  accurate  the  near-field  approach  will  actually  be  in 
solving  the  problem  of  source  radiation  in  this  complicated 
environment 


REFERENCES 

1.  Kinslcr,  Lawrence  E.,  Austin  R,  Frey,  Alan  B. 
Coppens,  and  James  V.  Sanders.  Fundamentals  of 
Acoustics.  New  York:  J.  Wiley,  2000,  Print. 

2.  Abraham,  D.  A.  "Array  Signal  Processing  for  Sonar: 
Short  Course."  Proc.  of  159th  Meeting  of  the 
Acoustical  Society  of  America,  Baltimore,  MD.  Print 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704-0788 

The  public  reporting  burden  for  this  colleclion  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructiorrs,  searching  exislmg  data  sources,  gathermg  and 
matnlainmg  the  data  rieedod.  and  compleling  and  reviewing  the  coflection  of  inlormation  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  the  burden,  io  the  Department  of  Defense,  Executive  Service  Directorate  (0704*0168)  Respondents  should  be  aware  thal  nolwilhslanding  any  other  provision  of  law,  no 
person  shall  bo  subject  to  any  penally  for  failing  lo  comply  with  a  collection  of  information  if  it  does  not  display  a  curranlty  valid  0MB  conirol  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION, 

1.  REPORT  DATE  pO-MM-yyYy^  2.  REPORT  TYPE 

4- 1 0-20 1 2  Final 

3.  DATES  COVERED  (Fmm  -  To) 

I  Feb  08 -29  Feb  I2 

4.  TITLE  AND  SUBTITLE 

Esiimation  of  Ship  Radiaicd  Noise  Based  on  Ncarficld  and  Structural  Acoustic 
Mcasurcmcnis 


6,  AUTHOR(S) 

Barber,  D.  Christopher;  Fowler.  Brian 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

NOOO  14.08. 1 -0427 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME{S)  AND  ADDRESS{ES) 

The  Pennsylvania  Slate  University,  Applied  Research  Laboratory,  P.  O.  Box  30,  State  College, 
PA  16804-0030 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


N/A 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Diane  Gales,  Code  ONR  BD025.  Office  of  Naval  Research,  875  North  Randolph  Street, 
Arlington,  VA  22203-1995 


10.  SPONSORyMONfTOR’S  ACRONYM(S) 

ONR 


11.  SPONSOR/MONtTOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 
Approved  for  Public  Release;  distribution  is  Unlimited. 


13.  SUPPLEMENTARY  NOTES 
N/A 


14.  ABSTRACT 

The  radiuted  noise  of  ships  or  other  underwater  sources  arc  typically  characterized  in  tenns  of  a  far-field,  plane-wave  equivalent  source  level  based 
on  measurements  assumed  to  have  been  acquired  in  a  free- field  environment  such  as  a  deep  water  test  range.  Measurement  of  ship  noise  in  a 
harbor  environment,  where  multiple  refiections,  high  background  noise  and  short  propagation  paths  are  the  nonn,  violates  the  conditions  that 
assume  the  ship  is  a  radiating  simple  source.  Careful  analysis  is  required  to  arrive  at  a  valid  estimate  of  far-field,  frce-field  source  levels  from  such 
measurements. 


15.  SUBJECT  TERMS 

Radiated  noise  measurements,  far-field.  wavefronts,  plane-wave,  harbor  environment,  propagation 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF 

PAGES 

5 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

D.  Christopher  Barber  /  David  L.  Bradley 

U 

U 

U 

uu 

19b.  TELEPHONE  NUMBER  (Inciudo  ama  code) 
814-863-9916 

Standard  Form  298  (Rev.  8/98) 

Presented  by  ANSI  Std.  239  18 
Adobe  Professtonel  7.0 


INSTRUCTIONS  FOR  COMPLETING  SF  298 


1.  REPORT  DATE.  Full  publication  date,  including 
day,  month,  if  available.  Must  cite  at  least  the  year  and 
be  Year  2000  compliant,  e.g.  30-06-1998;  xx-06-1998; 
xx-xx-1998. 

2.  REPORT  TYPE.  State  the  type  of  report,  such  as 
finai,  technical,  interim,  memorandum,  master's  thesis, 
progress,  quarterly,  research,  special,  group  study,  etc. 

3.  DATES  COVERED.  Indicate  the  time  during  which 
the  work  was  performed  and  the  report  was  written, 

e  g.,  Jun  1997  -  Jun  1998:  1-10  Jun  1996;  May  -  Nov 
1998;  Nov  1998. 

4.  TITLE.  Enter  title  and  subtitle  with  voiume  number 
and  part  number,  if  applicable.  On  classified 
documents,  enter  the  title  classification  in  parentheses. 

5a.  CONTRACT  NUMBER.  Enter  all  contract  numbers 
as  they  appear  in  the  report,  e.g.  F33615-86-C-5169. 

5b.  GRANT  NUMBER.  Enter  ail  grant  numbers  as 
they  appear  in  the  report,  e.g.  AFOSR-82-1234. 

5c.  PROGRAM  ELEMENT  NUMBER.  Enter  all 
program  element  numbers  as  they  appear  in  the  report, 
e.g.  61 101  A. 

5d.  PROJECT  NUMBER.  Enter  all  project  numbers  as 
they  appear  in  the  report,  e.g.  1F665702D1257;  ILIR. 

5e.  TASK  NUMBER.  Enter  all  task  numbers  as  they 
appear  in  the  report,  e.g.  05;  RF0330201;  T4112. 

5f.  WORK  UNIT  NUMBER.  Enter  all  work  unit 
numbers  as  they  appear  in  the  report,  e.g.  001; 
AFAPL30480105. 

6.  AUTHOR(S).  Enter  name(s)  of  person(s) 
responsible  for  writing  the  report,  performing  the 
research,  or  credited  with  the  content  of  the  report.  The 
form  of  entry  is  the  last  name,  first  name,  middle  initial, 
and  additional  qualifiers  separated  by  commas,  e.g. 
Smith,  Richard,  J,  Jr, 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND 
ADDRESS(ES).  Self-explanatory. 


8.  PERFORMING  ORGANIZATION  REPORT  NUMBER. 
Enter  all  unique  alphanumeric  report  numbers  assigned  by 
the  performing  organization,  e.g.  BRL-1234; 

AF  WL-TR-85-40 1 7-Vol-2 1  -PT-2. 

9.  SPONSORING/MONITORING  AGENCY  NAME(S) 
AND  ADDRESS(ES).  Enter  the  name  and  address  of  the 
organization(s)  financially  responsible  for  and  monitoring 
the  work. 

10.  SPONSOR/MONITOR'S  ACRONYM(S).  Enter,  if 
available,  e.g.  BRL,  ARDEC,  NADC. 

11.  SPONSOR/MONITOR'S  REPORT  NUMBER(S). 

Enter  report  number  as  assigned  by  the  sponsoring/ 
monitoring  agency,  if  available,  e.g.  BRL-TR-829;  -215. 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT.  Use 
agency-mandated  availability  statements  to  Indicate  the 
public  availability  or  distribution  limitations  of  the  report.  If 
additional  limitations/  restrictions  or  special  markings  are 
indicated,  follow  agency  authorization  procedures,  e.g. 
RD/FRD,  PROPIN,  ITAR,  etc.  Include  copyright 
information. 

13.  SUPPLEMENTARY  NOTES.  Enter  information  not 
included  elsewhere  such  as:  prepared  in  cooperation 
with;  translation  of;  report  supersedes;  old  edition  number, 
etc. 

14.  ABSTRACT.  A  brief  (approximately  200  words) 
factual  summary  of  the  most  significant  information. 

15.  SUBJECT  TERMS.  Key  words  or  phrases  identifying 
major  concepts  in  the  report. 

16.  SECURITY  CLASSIFICATION.  Enter  security 
classification  in  accordance  with  security  classification 
regulations,  e.g.  U,  C,  S,  etc.  If  this  form  contains 
classified  information,  stamp  classification  level  on  the  top 
and  bottom  of  this  page. 

17.  LIMITATION  OF  ABSTRACT.  This  block  must  be 
completed  to  assign  a  distribution  limitation  to  the  abstract. 
Enter  UU  (Unclassified  Unlimited)  or  SAR  (Same  as 
Report).  An  entry  in  this  block  is  necessary  if  the  abstract 
is  to  be  limited. 


standard  Form  298  Back  (Rev.  8/98) 


